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Abstract. We calculate dilepton spectra in heavy-ion collisions using a coarse-graining 
approach to the simulation of the created medium with the UrQMD transport model. This 
enables the use of dilepton-production rates evaluated in equilibrium quantum-field theory at 
finite temperatures and chemical potentials. 


Due to their penetrating nature, the invariant-mass and transverse-momentum spectra of 
dileptons are considered unique probes for the in-medium properties of hadrons in strongly 
interacting hot and dense matter, created in heavy-ion collisions [1, 2]. The dilepton-production 
rates are proportional to the electromagnetic current-correlation function of the source [3] and 
thus provide insight into the corresponding excitations carrying the quantum numbers of the 
electromagnetic current. Particularly from the restoration of chiral symmetry one expects 
substantial modifications of the spectral shape of this current-correlation function. Particularly 
one hopes to gain insight into the microscopic realization of the chiral-symmetry restoration, e.g., 
whether it is due to dropping in-medium hadron masses as predicted by the Brown-Rho-scaling 
conjecture or models implementing the vector manifestation of chiral symmetry or substantial 
broadening of the pertinent light-vector-meson spectral shapes due to hadronic interactions 
within the medium (“melting-resonance scenario”) as predicted by quantum-field-theoretical 
effective hadronic many-body theory (for a review see [1]). 

Since the dileptons are produced and emitted undisturbed by final-state interactions during 
the entire evolution of the hot and dense medium, first a good understanding of all relevant 
microscopic sources of electromagnetic radiation in this medium is important. 

For heavy-ion collisions at ultrarelativistic energies, in the early hot stage the relevant degrees 
of the medium are quarks and gluons forming a quark-gluon plasma (QGP) close to local thermal 
equilibrium, as is inferred from the successful description of the hadronic bulk observables like pT 
spectra and angular anisotropies Vn by hydrodynamical models for the evolving medium. Here 
the main source of dileptons is the quark-antiquark annihilation q + q^i^+i~ (with i = e, p). 
In the later stages of the fireball evolution (or, in the case of heavy-ion collisions at lower energies, 
the entire evolution) the medium is described as a hot hadronic gas. In the low invariant-mass 
range, 2me < M < IGeV the dilepton sources are thus hadronic and phenomenologically well 
described within a vector-meson dominance (VMD) model [4], i.e., the assumption that the 
electromagnetic current of the hadrons is proportional to the light-vector meson field. This also 
implies that the electromagnetic transition form factors for the various Dalitz decays of hadron 
resonances contributing substantially to the dilepton yield in the low-mass tail are describable 





Figure 1. The dimuon spectra in 158^4 GeV In+In collisions at the CERN SPS in (a) the 
low-mass (LMR) and (b) the intermediate-mass region (IMR) in comparison to the data from 
the NA60 Collaboration [5, 6, 7]. In the LMR the dilepton yield is dominated by the decays of 
the in-medium p mesons (dotted blue line) and in the IMR by QGP and multi-pion hadronic 
processes. The comparison with a dilepton rate based on a thermal-p contribution, neglecting the 
effect of baryon couplings, shows the importance of the these contributions for the tremendous 
broadening of the p-spectral function as well as the yield in the very-low-mass range close to the 
two-muon threshold. Microscopically these contributions originate from Dalitz decays of various 
baryon resonances. 


within the VMD model. In the intermediate-mass range also “multi-pion processes” become 
relevant. 

All these microscopic processes have to be evaluated in a hot and dense medium, which 
is at present available in detailed many-body quantum-field theoretical calculations only for a 
medium in thermal equilibrium, with some chemical off-equilibrium taken into account in terms 
of effective baryon densities and pion- as well as kaon-chemical potentials [8, 9]. It has been quite 
well established that the observed excess of low-mass dilepton production in heavy-ion collisions 
as compared to the expectations from the “hadronic cocktail” (i.e., from a simple scaling of 
the findings in pp collisions by the number of nucleon collisions in the nuclear reaction) can be 
understood by effective hadronic models that predict a tremendous broadening of the light-vector 
mesons’ spectral shape with small mass shifts [10, 11, 12, 13]. 

On the other hand, one also has to describe the evolution of the hot and dense medium 
in a reliable way. This is usually done either by employing macroscopic descriptions, reaching 
from quite simple thermal-fireball parametrizations (“blast-wave fits”) [14] to detailed ideal or 
viscous hydrodynamics simulations [15, 16, 17, 18, 12, 19], or microscopic transport simulations 
[20, 21, 22, 23, 24, 25, 26, 27, 28, 29]. 

In the latter approach it is challenging to fully implement the medium modifications of the 
partonic and hadronic dilepton sources. In our approach this problem is solved by mapping 
the kinetic description of the fireball evolution via Monte-Garlo transport simulations to a state 
close to local thermal equilibrium. In this work, an ensemble of collision events simulated by 
the UrQMD model is used to calculate a coarse-grained energy-momentum tensor and net- 
baryon-number current on a space-time grid. Then, using the Eckart definition of the local rest 
frame, a temperature and baryo-chemical potential can be extracted via an equation of state 
(EoS), making use of a description of the energy-momentum tensor in terms of an anisotropic- 
fluid ansatz [30], where the pressure anisotropy, particularly relevant in the early stages of the 
fireball evolution, is taken into account. For the EoS we match a parametrization of lattice-QGD 
results [31] with a hadron-resonance gas EoS [32], which match very well in the region of the 









Figure 2. (Color online) (a) Dielectron invariant-mass spectrum for Ar-|-KCl collisions at 
Eiah = 1.76AGeV and (b) for Au-|-Au at Fiab = 1.23AGeV. The results are normalized to the 
average total number of per event and shown within the HADES acceptance. The results 
for Ar-|-KCl are compared to the data from the HADES Collaboration [35]. The comparison 
of the dilepton spectra between using the full model and omitting the contributions of direct 
baryon couplings to the p and oj self-energies underlines the prevalence of the corresponding 
baryon-resonance Dalitz-decay processes processes for the in-medium modifications of these 
vector mesons, shifting substantial spectral strength to the low-invariant-mass region. 


parton-hadron transition around the critical temperature of — 170 MeV. 

In recent work [33, 34], we have combined this coarse-graining method for the bulk-medium 
evolution of the strongly interacting fluid via the UrQMD transport model with state-of-the art 
in-medium dilepton-production rates [8, 9] in order to evaluate the dilepton invariant-mass and 
Pt spectra in 158AGeV In-|-In collisions as measured by the NA60 Collaboration at the CERN 
SPS and in 1.76AGeV Ar-|-KCl and 1.23AGeV Au-I-Au collisions as measured by the HADES 
collaboration at the GSI SIS. 

As exemplified by the comparison of the so calculated dimuon invariant-mass spectra in 
158A GeV In-I-In Collisions with the NA60 data (Eig. 1) and corresponding dielectron invariant- 
mass spectra for 1.76AGeV Ar-|-KCl collisions with data from the HADES collaboration (Eig. 2, 
left panel), the model successfully describes the data. In the right panel of Eig. 2 we additionally 
provide our prediction for the 1.23AGeV Au-|-Au collisions as also measured by the HADES 
collaboration. 

Both figures underline the importance of the interaction of the light vector mesons with 
baryons in the medium [in the present model [1] N(938), A(1232), N(1440), N(1520), A(1620), 
A(1700), N(1720), A(1900), N(2000)]. As expected, these processes are prevalent at the low 
GSI-SIS energy with high net-baryon densities but also relevant at CERN-SPS energies and even 
at higher RHIC energies. This is due to the fact that the strong interaction is CP invariant and 
thus not the net-baryon density, Nb — N-^, but the total baryon-antibaryon density, Nb + 
is relevant for the in-medium modifications of the light-vector mesons’ spectral distributions. 

In the near future we also plan to apply the model to dilepton production at RHIC energies 
as well as photon production in heavy-ion collisions. 
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